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ABSTRACT 

There is a large observational scatter toward low velocities in the stellar mass Tully- 
Fisher relation if disturbed and compact objects are included. However, this scatter 
can be eliminated if one replaces rotation velocity with So. 5, a quantity that includes a 
velocity dispersion term added in quadrature with the rotation velocity. In this work 
we use a large suite of hydrodynamic N-body galaxy merger simulations to explore a 
possible mechanism for creating the observed relations. Using mock observations of the 
simulations, we test for the presence of observational effects and explore the relationship 
between So. 5 and intrinsic properties of the galaxies. We find that galaxy mergers can 
explain the scatter in the TF as well as the tight So.5-stellar mass relation. Furthermore, 
So. 5 is correlated with the total central mass of a galaxy, including contributions due to 
dark matter. 

Subject headings: galaxies: formation — galaxies: interactions — galaxies: kinematics 
and dynamics — galaxies: evolution 



1. Introduction 



One of the goals of a theory of galaxy formation is to reproduce the scaling relations of both 
early and late type galaxies. There has been significant work on understanding galaxy formation 
from ana lytic, semi-analytic , and numerical perspectives dating back to the collapsing gas cloud 
model of lEggen et al.l (|l962l ). For early type galaxies, scaling relations between their sizes, lumi- 
nosities or stellar masses, and velocity dispersions form a two-d imensional plane called the "Fun- 



damental Plane" (Dressier et al 



1987 



Diorgovski &: David 119871 ). For late type galaxies, a simple 



plane does not exist (jBurstein et al.lll997l ). and the relatioi i between luminositie s or stellar masses 
and rotation velocities (i.e., the Tully-Fisher (TF) relat ion; iTully Sz Fisheiill977l ) is independent of 
size and surface brightness (e.g., ICourteau &: Rixlll999l ). 



The Fundamental Plane was found to be a consequence of the virial theorem, if one assumes a 
constant mass-to- light ratio and homologou s mass structure (jFaber et al.lll987l : lDiorgovski &: Davis 
19871 : iBender et al.lll993l : lPahre et al.lll998l ). However, the situation is more complicated since there 
is a tilt observed between the Fundamental Plane and the plane resulting from the virial theorem. 
This tilt is likely due to a combination of s t ellar population effects and another factor which has 
yet to be determined ( e.g., 



Pahre et al 



1998 



Gerhard et al. 



2001 



Bertin et al. 



2002 



Hopkins et al 



20081 ). iHopkins et al.l ( 20081 ) argue that this other factor results from a varying ratio of total enclosed 



mass within the effective radius (Re) to stellar mass. This ratio must depend on mass such that 
lower mass ellipticals are more baryon-dominated within than their higher mass counterparts. 



For disk galaxies, the scaling relations between velocity, luminosity or stellar mass, and size can 
broadly be understood by semi-analytic galaxy formation models which are based upon the dissi- 
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pationless collapse of a cold dark matter halo and angular momentum conservation (e.g., 


Mo et al. 


1998; 


van den Bosch 


2001; 


Dalcanton et al. 


1997: 


de Jone & Lacev 


2000 


; Button et al. 


20071). Such 


mode 


s have their roots in the analytic formulatio 
) which in turn have their roots in the wo 


ns of 


White & Prenk ( 1 


991 


) and 


Fall & Efstathiou 


(1980 


rks of 


Silk 


( 


1977 


) and 


Rees & Ostriker ( 


1977). A 



detailed prediction of disk galaxy scaling relations, which necessitates simultaneous matching of the 
zero-points o f the TF relation, galaxy sizes, and the luminosity /stellar mass functions has proven 



elusive (e.g., iBenson et al 



2003: 



these properties (IDutton Sz van den Bosch 



Dutton et al.ll2007l ). Models can be constructed that reproduce 
2009! ). but they req uire efficient mass outflo ws (a.k.a. 



feedback), and the absence of dark halo adiabatic contraction (jBlumenthal et al.l 119861 ). On the 
purely numerical front, hydrodynamic + N-body codes have been able to form disk-like systems 
only under controlled conditions and for halos carefully s e lected based upon their merger histories 



(e.g. , iGovernato et al.l 12004 120071 : [Robertson et al.l |200J; IScannapieco et al.l l2009l ) . Furthermore 



the format ion of realistic bulg e-less galaxies remains an unsolved problem for cosmological simula- 
tions (e.g. jMaver et ahlbood ). 



Cosmological hydro dynamical simulations have been able to form galaxies with disks that 
reproduce the slope of the TF relation (e.g. Navarro & Steinmetz 2000; Eke et al. 2001, Portinari 
et al. 2007, Governato et al. 2007; Piontek & Steinmetz 2009). However, as with analytic and 
semi-analytic models, reproducing the TF zero point, and observed disk sizes has been a challenge. 
If the rotation velocity is measured at a radius which encloses most of the baryons, then simulated 
galaxies can reproduce the zero point of the TF relation (e.g. Governato et al. 2007). However, 
if the rotation velocity is measured at a smaller radius (such as 2.2 K-band disk scale lengths), 
the same simulations of Governato et al. (2007) over-predict the rotation velocities (Dutton & 
Courteau 2008). This failure is due to the simulations having the wrong distribution of baryons 
and/or dark matter within the optical parts of galaxies. 

In summary, although the scaling relations of early and late type galaxies can be explained in 
broad terms, detailed predictions have yet to fully mature. In addition, since early type galaxies 
are generally expected to form through mergers of late type systems, the scaling relations of early 
types should therefore descend from those of late ty pes. There has been work with simulations of 



mergers of late type s to try to und e rstan d this (e.g.. iRobertsoii et al.l l2006l : ICovington et al.l 120081 : 



Hopkins et ahlbood ). IPekel fc Co^] (j20od ). iHopkins et alJ kood ). and Covington et al. (in prep) in 



particular argue that late type scaling relations may evolve into those of early types if the varying 
gas content of disk galaxies with mass and dissipation during mergers are taken into account in 
models. However, not all galaxies in today's universe are the well-ordered disk or elliptical systems 
for which most theories make predictions. Many galaxies have disturbed or irregular morphologies 
and /or kinematics, and the fraction of these galax ies increases with redshift to at least z ~ 1 (see 



e.g.. Ivan den Bergh et al.lll996l ; iKassin et al.ll2007l. for morp h ologie s and kinematics, respectively). 
With this in mind, we examine the findings of iKassin et al.l (j2007l ) which attempt to include such 
irregular systems in a study of the Tully-Fisher scaling relation, and to link early and late types 
through their kinematics. 
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In particular , we focus on the sca ling laws that relate galactic orbital speeds to stellar masses: 
the TF relation (jTullv &: Fished 119771 ) for rotating disk gal axies, which relates stel lar mass to ro- 
tation velocity (Vrot); and the Faber- Jackson relation (FJ) (jFaber &: Jacksonlll976l ) for spheroidal 
(elliptical) galaxies, which relates stellar mass to velocity dispersion (o"). T he conventionally as- 
sume d pathway for forming ellipticals is through mergers of disk galaxies (iBarnes &: Hernquist 
1992), in which case the FJ relation descended from the TF relation, as modified by mergers. 



The TF and FJ relations relate the specific kinetic energy required to support a galaxy to the 
attractive potential well create d by the mas s of th e galaxy. Using the DEIMOS multi-object spec- 
trograph at Keck Observatory, iKassin et al.l (|2007l ) observed the kinematics and masses of emission 
line galaxies in a redshift range from 0.1 to 1.2 from the All Wavelength Extended Groth Strip 
International Survey (AEGIS), incorporating morphologically normal and disturbed galaxies. They 
find that the disturbed galaxies do not lie on the standard TF relation which only considers global 
rotational support. Many of the disturbed galaxies have rotation velocities that are well below the 
normal values. This morphological dependence of the TF scatter is well known. St udies of local 



galax ies have found larger scatter for close pair s of galaxies with kinern atic disorders (IBarton et al 



200lh and when peculiar galaxies are added ([Kannappan et al.ll2002l ). At z ~ 0.5, iFlores et al 
(j2006l ) find larger TF scatter for galaxies wit h "complex k inema tics," which was confirmed by the 
larger sample and more detailed analysis of iPuech et al.l (|2008l ). This correlation of scatter with 
morphology is not surprising, as one might expect kinematic scaling relations based on pure rota- 
tional support to fail for cases where the systems are not in globally ordered motion. Galaxies that 
are merging, disturbed, or rapidly evolving may exhibit disordered kinematics, so that considering 
only rotation will not capture all the kinematic support that balances the gravitational attraction 
of the potential. 



Interestingly, when lWeiner et al.l (|2006l ) and lKassin et al.l (|2007l ) incorporated the random in- 
ternal motions of galaxies, as measured by velocity dispersion, along with rotational velocity into a 
new parameter, Sq.s = ^0.5V^^^ + 0"^, then all galaxies regardless of morphology fell onto a single 
relation between So. 5 and stellar mass. This new parameter allowed for all galaxies to be incor- 
porated in a measurement of the TF relation, instead of requiring the TF to be restricted to a 
sub-sample of galaxies selected to be morphologically and/or kinematically well-or dered. In add! 



tion, t his new TF relation is consistent with the FJ relation for elliptical galaxies from lGallazzi et al 



(|200d ). It is perhaps surprising that a set of non-homologous systems fall on the same kinematic 
relation, suggesting that the relation contains information on the kinematic evolution of galaxies. 

Given the complexity of effects at work, the constraints that produce this tight kinematic re- 



lation are r iot clear. The veloc ity dispersion as measured in gas-phase kinematics by IWeiner et al 



(120061 ) and iKassin et al.l (120071 ) appears to provide kinematic support similar to a thermalized ve- 
locity dispersion, but it probably represents spatially unresolved velocity gradients rather than 
truly randomized motions. Observational effects, such as blurring, can transform chaotic motions 
or spatially unresolved rotation into apparent dispersion. The So. 5 quantity is formulated to mea- 
sure a total specific kinetic energy; although these systems are non-equilibrium, it is likely that the 
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virial theorem remains a good approximation for the balance between kinematics and gravitational 
potential. In order to investigate these possibilities and to explore the origin of the So.5-stellar mass 
relation, we analyze a large suite of hydrodynamical galaxy merger simulations. Mock observations 
of this simulation set provide kinematic data on isolated disk galaxies, disturbed galaxies, merg- 
ing/overlapping galaxies, and elliptical merger remnants. The simulations allow for a comparison 
with observational results as well as a parallel analysis of intrinsic kinematics and observational 
effects. 



Methods 



2.1. Galaxy Merger Simulations 

We exploit a large suite of binary galaxy merger simulations originally used to study feedback 
and star-formation in galaxy mergers. It provides a rich source of kinematic data for galaxies with 
a wide range of morphologies. For this study use a selection of 50 merger simulations from this 



suite . For full details on the merger si mulations we refer the reader to previous studies (|Cox et al 



2004 



Cox 



2004 



Cox et al. 



2006 



20081 ). but we include a brief description here for completeness. 



The numerical simulations performed in this work use the N-Body/SPH code GADGET 
([Springel et al.ll200lh . Hydrodynamics are included via the Lagrangian techni que of smoothed parti- 
cle hy drodynamics (SPH). We use the "conservative entropy" version of SPH (jSpringel &: Hernquist 
2002). Gas is assumed, for simplicity, to be a primordial plasma that can radiatively cool via atomic 
and free-free emission. 

All of the numerical simulations presented here include star formation. Stars are formed in 
regions of gas that are above a critical density for star formation at a rate proportional to the 
local gas density and inversely proportional to the local dynamical time-scale. The efficiency of 
star formation is fixed by requiring star f ormation to follo w the observed correlation between gas 
surface densities and star- formation rate (|Kennicuttlll998l ). 



We also include a simple prescription to simulate the effects of feedback from massive stars. 
This feedback acts to pressurize the interstellar medium and reg ulates the c o nvers ion of gas to stars. 
Details of this model and the parameter choices can be found in lCox et al. ( 2006 ). Specifically, the 
simulations studied in this paper used the n2med parameter set. Under these assumptions the gas 
pressure increases as the density squared; i.e, star-forming gas has a "stiff" equation of state. 

The simulations adopt a gravitational softening length of /i = 400 pc for the dark matter 
particles and 100 pc for the stellar and gas particles. We remind the reader that, in GADGET, forces 
between neighboring particles become non-Newtonian for separations < 2.3 times the gravitational 
softening length. 



The disk galaxy models and orbits are cosmologically motivated, but the simulations are not 
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cosmological since the two galaxies are isolated. All mergers are binary mergers of disk galaxies. 
Rather than being produced by previous mergers, the progenitor disk galaxies are constructed 
in equilibrium and contain dark matter, an exponential stellar disk, an extended exponential gas 
disk, and some contain a dense central bulge. The dark matter halos are assumed to have an NFW 
profile, whereas the stellar bulges have a three-dimensional exponential profile. Bulge to disk stellar 
mass ratios vary from 0.02 to 0.25. The portion of the suite used in this study contains two types 
of progenitor galaxy models: 



1. "Sbc" galaxies, which are modeled after local Sbc-type spirals, with a small bulge and an 
initial gas mass roughly equal to the initial stellar mass. 

2. "G" galaxies, which span a range of mass, bulge fraction, and gas fraction. Their properties 
are take n from statistica l samples of local late type galaxies, including the Sloan Digital Sky 
Survey dVork et al.1 boool l . 



The Sbc mergers are all major mergers of identical progenitors with a stellar disk of 3.92 x 
1O^°M0, a gas disk of 5.36 x IO^^Mq, a stellar bulge of 1.00 x 10^°Mq, and a dark matter halo of 
81.2 X lO^'^M©. The Sbc series contains the widest variation of initial merger orbits and orienta- 
tions, with 17 different orbits run. This includes orbits with a wide range of pericentric distances, 
eccentricities between 0.6 and 1.0, and a variety of progenitor orientations, including prograde- 
prograde, prograde-retrograde, retrograde-retrograde, and some cases inbetween. The G series 
mergers include both major mergers between identical galaxies an d mergers of galaxies wi t h mas s 
rati os between 1:1 and 1:5 0. For more detail on these models see ICoxl ()2004l ). I Cox et al.l ()2008l ). 
and ICovington et al.l (j2008l ) . 



2.2. Mock Observations of the Simulations 



In order to make comparis ons with the observa tional results, we must first mock 'observe' the 
galaxies in the sarn e manner as iKassin et al.l (j2007l ). They used an algorithm called ROTCURVE 
(jWeiner et al.l 120061 ). This algorithm was designed to obtain kinematic measurements for as many 
galaxies as possible by allowing model seeing-compensated rotation curves to be fit using only 
the 2-d spectra. ROTCURVE creates an intrinsic model for the emission intensity, velocity, and 
dispersion, blurs this model to simulate seeing, and fits the model to the data. The intrinsic model 
used is: 

{v-V{x)f. 



I{x, v) = G{x) exp(- 



G{x) = 
V{x) 



hot , {x 
exp( 



27rri 



2r} 



Vrox.— arctan(x/rv), 
vr 



(1) 
(2) 
(3) 
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where G{x) is an assumed Gaussian light distribution along the slit, V{x) is the rotation curve, with 
asymptotic velocity Vrot and knee radius r^, and cJ2d is the velocity dispersion, which is assumed to 
be constant along the slit. In order to fit the spectra, v) is blurred using a 1-d Gaussian in the 
spatial direction. After blurring, moments are taken of the blurred model in the velocity direction 
in order get values of average velocity and dispersion for each bin within the slit. These values are 
then compared against the observed rotation curve by computing a value, is then minimized 
via adjustments of Vrot and a2d- The other potentially adjustable parameter, r^, is held fixed at 0.2 
arcsecs. Observational results do not change significantly if is varied within a reasonable range 
of 0.1 to 0.3 arcsec. 

In order to analyze the simulations, we use an analogous algorithm. Each simulation has snap- 
shots every 40-100 Myr, and each of these snapshots is viewed from 11 different evenly distributed 
angles. The simulations are assumed to be at three different redshifts {z ~ 0.03, z ~ 0.3, and 
z ~ 1-0), which approximately cover the redshift range in Kassin et al. We include the lowest 
redshift bin to examine the relation with very little blurring. This results in a sample of roughly 
150,000 mock observations. Physical distances are converted to arcsecs given the assumed redshift 
of the observation, using cosmological parameters h = 0.7, Qm = 0.3, and = 0.7. The effects of 
seeing are taken into account by smearing particle positions using a Gaussian with a = 0.3 arcsecs. 
If the projected separation of the two galaxies is less than 0.7, arcsecs then the two galaxies are 
treated as one object. This allows us to include merging and overlapping galaxies in the sample. 
We orient the mock slit by fitting an ellipse to the simulated stellar density profile. All mock 
observations use a slit width of 1.0 arcsec as in the observations. The mock slit length varies with 
redshift beginning with a length of 4.0 arcsecs at z ~ 1 and scaling with physical galaxy size, as in 
the observations. 

Within the slit, stellar particles - either all stars or just new stars (see Section 14. 2p - are 
separated into 0.1 arcsec bins, mimicking the pixel width in the observations. Rotation curves are 
calculated using the average velocities and dispersions of the particles in each bin. We fit G( x) to 



the bin stellar mas ses and calculate I{x, v) and V{x) as above. Unlike IWeiner et al.l (|2006l ) and 



Kassin et al.l (|2007l ). we do not hold fixed, but instead assume a value of 0.2 arcsecs at z ~ 1 and 
scale that value with the physical size of the galaxy. If we do not scale and the slit length with 
redshift, we get spurious systematic offsets in Vrot when comparing the lowest redshift bin with the 
two higher redshift bins. Inclination is also calculated, using the thin disk approximation as in the 
observational work, and is used to correct Vfot to edge-on. We remove all galaxies from the sample 
with inclinations less than 50 degrees, as we find that the inclination fits become unreliable below 
that value. Similarly, Kassin et al. remove the galaxies with inclinations less than 30 degrees. 
The observational sample is also cut for nearly edge-on disks (inclination greater than 70 degrees) 
because of dust effects on stellar mass estimates. Since the simulations do not suffer from this same 
limitation, we do not cut at large inclinations. 
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3. Comparison of the Simulations with Observations 

3.1. Evolution of a Single Galaxy Merger 

To illustrate how the kinematics of galaxies change during a merger event, we follow the time 
evolution of a single galaxy merger. Even though this is only a specific example case from the set 
of simulations studied, it is typical of the mock observations of the other simulations, as evidenced 
by the results of Section 13. 2i The simulation used for this study is a merger of two identical Sbc 
galaxies. Initially, the galaxies are set on a parabolic prograde-prograde orbit with a pericentric 
distance of 11.0 kpc. We depict the evolution of kinematic quantities during the merger viewed 
from a single angle (Figure 1), and follow the galaxies' evolution on the kinematic relations (Figure 
2, numbers). The merger is assumed to be at a redshift of one. 

We start off at time=0 Gyr with the progenitors, which are designed to lie on the observed TF 
ridge-line. The progenitor is demarcated as '1' in Figures 1 and 2. We follow one of the galaxies 
until the projected distance between the galaxies becomes less than 0.7 arcsecs after which we 
follow the merger as a single system. The merging galaxies have two encounters and merge during 
the second encounter. The first encounter occurs at 0.6 Gyr after the start of the simulation. 
This encounter results in a significant increase in the velocity dispersion by ~ 60 km s^^ and a 
sharp decrease in the rotation velocity by ~ 55 km s~^. During this encounter, the galaxy appears 
morphologically disturbed, especially in its outer disk. Shortly following the encounter, at 0.7 Gyr, 
the galaxy is a slightly low rotation velocity outlier in the TF ('2' in Figures 1 and 2). Between 
the first and second/final encounters, the rotation velocity of the galaxy gradually increases by 
~ 20 km s~^ as tidal debris that was removed from the disk settles back. The final coalescence 
occurs at 1.75 Gyr after the start of the simulation ('3' in Figures 1 and 2). At coalescence, the 
galaxy's rotation velocity drops abruptly to a meager 20 km s~^. The mock observational slit now 
contains both the first progenitor and the remains of the second progenitor, which has been severely 
disrupted. Projection results in the overlapping of the two progenitors within the slit and pollutes 
the rotation curve, causing the dramatic decrease in rotational velocity. Consequently, the system 
lies furthest from the TF law at this time ('3' in Figures 1 and 2). At the final coalescence, we 
find an increase in the velocity dispersion by ~ 100 km s~^. From this point on, the kinematics are 
dominated by dispersion rather than rotation. The final stage ('4' in Figures 1 and 2) depicts the 
merger remnant, which is a rotating elliptical galaxy. In contrast to the other kinematic quantities. 
So. 5 remains relatively stable throughout the merger process, increasing primarily when the two 
progenitors combine to create a more massive system. 

The tendency for rotation velocity to decrease as velocity dispersion increases during galaxy 
interactions and mergers provides a simple mechanism to explain the observed So.s-stellar mass 
relation and large scatter in the conventional TF. We note, however, that observed So. 5 is not 
strictly conserved, even when mass is constant. During a close encounter all of the 'lost' rotation is 
not immediately converted into velocity dispersion. Rather, the tendency during such an encounter 
is for the apparent rotation velocity to briefly drop to very low values while the velocity dispersion 
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increases only modestly. The effects of this lack of strict conservation of S0.5 fall within the scatter 
of the So.5-stellar mass relation. 



3.2. Kinematic Relations for an Ensemble of Simulated Galaxy Mergers 

In order to further test this explanation for the observed kinematic relations, we plot rotation 
velocity and So. 5 versus stellar mass for 500 randomly chosen s napshots from set o f 50 merger 



simulations and compare to similar plots for 544 galaxies from iKassin et al.l (|2007l ) (Figure 2). 
Each of these snapshots is viewed from an angle chosen at random from the 11 different evenly 
distributed viewing angles. These 500 snapshots are analyzed assuming z ~ 0.03, z ~ 0.3, and 
z ^ 1.0. We find that at all redshifts the simulated TF has a significant scatter to low T4ot 
(spanning ~ 1.5 dex), similar to the observations. Furthermore, this scatter is correlated with close 
encounters. Sixty-four percent of the cases with scatter greater than 0.5 dex are encounters /mergers, 
whereas encounters/mergers are only 30% of the entire sample. The primary difference between 
the simulated and observed results is in the correlation between scatter and redshift. For the 
observations, scatter in the TF significantly increases at high redshift, whereas for the simulations, 
the scatter in the TF does not correlate strongly with redshift. This suggests that the observed 
evolution in scatter of the TF represents a real change in kinematics with redshift as opposed to 
an observational blurring effect. 

If we incorporate random internal motions using S0.5, then we find a single, relatively tight 
relation with stellar mass. Comparing different redshift bins we see that the slope fitted to the 
simulated galaxies stays within a range of ~ 2.8 — 2.9, whereas the observed slopes are typically 
slightly steeper but range between ~ 2.4 — 3.3. There is no significant trend in slope with redshift 
in either the observations or simulations. In the two higher redshift bins the scatter in So. 5 for the 
simulations is 0.08-0.09 dex. The intrinsic scatter in the observed relation ranges between 0.08 and 
0.12 dex and also has no systematic trend with redshift. The lowest redshift bin for the simulated 
results does show more scatter (0.16 dex), but this is likely due to difficulty in fitting unblurred 
rotation curves of disturbed objects, which often do not look like th e idealized ar c tan r otation 



curve. This is also outside of the range of the observations discussed in lKassin et al.l (j2007l ). 



While a detailed quantitative comparison of the samples is not warranted because the galaxies 
in the simulation suite are not a statistically representative sample of the types of galaxies that one 
would expect to find in the real Universe, the agreement with the observed relations demonstrates 
that galaxy mergers and interactions are a mechanism that can create large scatter in the TF 
concurrent with low scatter in the So.s-stellar mass relation. 

Furthermore, the simulated results show that the merging process moves galaxies up the So.s- 
stellar mass relation toward higher values of So. 5, and therefore that the relation ties together 
progenitors (spirals), merging galaxies, and merger remnants (ellipticals). The robustness of this 
relation has two causes. First, throughout the merger process there is a tendency for rotation and 



- 10 - 



velocity dispersion (or, random internal motions) to be anti-correlated. It has been known for a 
long time that galaxy mergers a re an effective mechanisn i for converting ordered rotational support 
into random pressure support (jToomre &: Toomrelll972l ). As noted above, it is surprising that a 
kinematic relation would hold for even merging cases. This results from a subtly different cause. 
Much of the apparent rotation velocity may be lost in an encounter when the two galaxies overlap 
in projection. This is because the two galaxies' orbital velocities and internal velocities rarely add 
coherently when randomly overlapped. However, velocity dispersions are not so easily lost. In fact, 
overlapping two galaxies with interfering rotation curves will increase velocity dispersion. Thus, 
when such an overlap occurs, apparent rotation typically decreases, whereas velocity dispersion 
typically increases. 



4. Observational Effects on the So.5-Stellar Mass Relation. 

From the observations alone, it is unclear the extent to which the So.5-stellar mass relation is a 
result of observational effects or is telling us something about the intrinsic properties of the galaxies. 
In this section, we explore the various observational effects on the kinematic measurements of the 
simulations and determine whether they can contribute to the observed relation. 



4.1. Blurring 

As galaxies are observed at higher redshifts, their angular extent becomes smaller. Conse- 
quently, the typical ~ 0.7" seeing produced by the atmosphere in the Kassin et al. observations 
effectively moves the stars from one spectral bin into other nearby bins within t he slit. This reduce s 



the apparent rotation velocity and produces an artificial velocity dispersion (jWeiner et al.ll2006l ). 
The fitting algorithm used by ROTCURVE takes this into account by fitting a blurred model and 
extrapolating back to an intrinsic model. However, an interesting question is how successful this 
procedure is, and whether or not blurring conserves S0.5. If S0.5 is conserved by blurring, then it is 
possible that blurring is partially responsible for the scatter in the TF and the tight So. 5 relation. 



In order to explore the effects of blurring, we take the fiducial merger simulation from ^3.11 
(the Sbc) and analyze the kinematics after blurring the particles using gaussians with cr's of 0.0", 
0.1", 0.3", 0.6", 0.8", and 1.6", assuming the galaxy is at a redshift z ~ 1 (Figure [3|). We use the 
four example snapshots from Figure 1: the undisturbed disk, disturbed disk, merger, and merger 
remnant. The effective blurring within the AEGIS sample is all between 0.1" and 0.3", given seeing 
of 0.7" (defined as full- width at half maximum). There is little difference between the unblurred 
kinematics and 0.1" kinematics for all snapshots shown. However, for both the disks and the 
remnant there is a decrease in Vfot and increase in a by the time we have reached 0.3". However, 
for all cases, the decrease in T4ot < 30%. Also, for all the blurring becomes much greater 

than i?5o, the radius that contains 50% of the stellar mass, the ROTCURVE algorithm breaks down 
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and Vrot drops to very low values while a rises. 

For the undisturbed disk, disturbed disk, and merger, S0.5 is approximately conserved by 
blurring as T4ot and a adjust in lockstep. However, for the remnant, by 0.3" blurring, S0.5 has 
already increased by ~ 50%. This results from a steeply-peaked central velocity dispersion. For 
zero blurring, the assumption of a constant a underestimates the kinetic energy contribution from 
random motion because the central peak is small with respect to the size of the slit. As this peak 
gets blurred, the fit improves. 

Therefore, it is true that as blurring becomes very large (>> R50), the ROTCURVE fitting 
procedure will break down and the rotational component of the kinematics will transform into 
dispersion. Also, for three of the four cases shown, S0.5 is conserved throughout this process. 
Therefore it is true that blurring could cause some additional spread in the TF and that using S0.5 
rather than rotation velocity would help to reduce this spread. However, for the range of blurring 
and galaxy sizes in the observational study, the ROTCURVE model performs relatively well, and 
certainly does not produce large enough errors to low 14ot to account for the observed scatter. This 
is also confirmed by the apparent lack of systematic changes in V^ot with redshift in Figure 2. Thus 
blurring cannot be the main reason for the observational findings. 

4.2. Tracking of Emission Lines 

A potential concern with our analysis is that the AEGIS observational results rely on Oil 3727 
emission line spectra from hot gas, whereas the kinematic analysis of the simulations described 
above uses the star particles. Specifically, the line emission is known to be correlated with areas 
of new star formation. Thus, we repeat the analysis from ^3.21 using only new stars that form in 
the simulation and assuming z = 1 (Figure S]). This gives us many fewer particles to work with, 
especially in the early stages of the merger when few new stars have formed. In order to provide 
sufficient statistics within the slit, we restrict the analysis to snapshots that have greater than 500 
new star particles. This removes some of the sample from early on in the simulations, but will 
provide a check for whether or not limiting the analysis to new stars affects the results. 

The TF and So. 5 relations are qualitatively quite similar to those calculated using all star 
particles, showing large scatter for TF and relatively small scatter for So. 5. This suggests that 
the observational difference between tracing kinematics using hot gas and stars does not have a 
significant effect on the results. There is some additional scatter in the new star So. 5 relation (0.14 
dex as compared to 0.08 dex). However, one must also remember that dust, which is not taken into 
account in the simulations, creates a counteracting effect. Within the simulations, analyzing only 
new stars will emphasize the central portions of the galaxies (where most stars form), however, 
these also tend to be the portions of the galaxy most enshrouded by dust. Thus, the real difference 
between observing emission lines (gas) or absorption lines (stars) is likely to be less. In order to 
truly pin down the effect of dust, one would have to create artificial spectral kinematics by fully 
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modeling the radiative transfer through the gas and dust. This wiU soon be possible using a newer 
version of SUNRISE (Jonsson, 2006 and in preparation). 



5. Intrinsic Kinematic Quantities 



We have explored some of the observational effects on the kinematic relations, however, we 
can further utilize the simulations to examine intrinsic kinematic quantities in order to better 
understand the origin of the So.s-stellar mass relation. If one assumes a tracer population within 
an isothermal sphere at a distance r >> rcore; where r^ore is the core radius of the isothermal sphere, 
and then allows that population to assume complete rotational support or complete support from 
random motion then one finds that Vdrc ~ V^cr, w here VcWc is the velocity with pure circu lar support 



and a is the dispersion for pure random support (jWeiner et al. 



2006 



Kassin et al.ll2007l ). This was 



the motivation for the multiplying factor of 1/2 in front of 14ot in Sq.s- One might expect that 
a smooth transition between rotational and random support, while conserving the density profile, 
would result in conservation of So. 5- Therefore, an important question is whether this idealized 
picture is correct and Sq.s really does trace the mass distribution. If so, this would provide a means 
of estimating total masses of galaxies including the contribution from dark matter. 



5.1. Kinematics as a Function of Radius 

We begin our analysis of intrinsic galaxy properties by examining kinematics as a function 
of radius. The So. 5 parameter is motivated by results obtained from the Jeans Equations for an 
isothermal sphere. However, our systems are not isothermal spheres, so we take a step back to see 
where the assumptions could be wrong. For steady state spherically symmetric systems with no 
rotation, the following relation holds (Binney and Tremaine, 1987, eq. (4-55)): 

r \dlnrdlnr I 

where Cr is the radial velocity dispersion, u is the mass density, P = 1 — a'^/a'^ is the velocity 
anisotropy parameter, and erg and cr^ are the velocity dispersions in the 6 and (f> directions. At a 
given radius, this can be represented as the simpler 

KL = kal (5) 

For the singular isothermal sphere, the second two terms in the parentheses are zero and we 
are left only with V^^j.^ = —a'^{d In u/d In r). For an isothermal sphere with a core radius rcore 
and r >> rcore the density term (d In u)/{d In r) ~ —2, giving us the constant. A; ~ 2, in the S0.5 
relation. For our simulation set, and for real galaxies, we are neither guaranteed that the density 
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term will be ~ —2 nor that the other two terms will vanish. Thus we look in detail at the values 
these terms are likely to have. 

The velocity dispersion is generally a weak function of radius. Thus the dispersion term, 
d{ln (T^)/d(ln r) is likely to be relatively small. The drop in dispersion with radius is quite small 
in the progenitors and is largest for the remn a,nts, where s i gnific ant dissipation and star formation 



may have steepened the profiles. As shown in lDekel et al.l (|2005l ) the dispersion in the remnants is 
relatively well represented by o" oc r""'^. Thus the contribution of the dispersion term in equation [J] 
is likely to be —0.2 or smaller, which is quite small compared to the presumed value of the density 
term. Because of tidal disturbance and dissipation, the stellar orbits in the outer portions o f 



the merger remnants do exhibit relatively large anisotropics with /3 ~ 0.2 — 0.5 (jPekel et al.ll2005l ). 
Furthermore this term is multiplied by a factor of two bringing the overall contribution to ~ 0.4—1.0. 
For progenitors, which typically have low anisotropies, one would expect the contributions from 
these terms to roughly cancel, whereas for remnants there would be a net contribution to k of 
~ 0.2 - 0.8. 

The largest of the three terms is the density term. This term has a value of roughly two for 
an isothermal sphere and three or four for the outer portion of an NFW or a Hernquist profile. 
Thus we very well might expect a steeper density relation, p oc r~", and consequently higher value 
of k than is found for the isothermal case. For each simulation we fit a value of a to the total 
mass density curve in the vicinity of R50, the radius that contains 50% of the stellar mass. For 
the simulations, the value of o at R50 ranges between ~ 1.5 — 3.5, with a peak around negative 
two (see Figured]). The slope is typically steeper for massive progenitors and for merger remnants. 
In summary, for progenitors, the density term is likely to be ~ 2, while the other terms will be 
small. For remnants, the density term will often be closer to 3, but the anisotropy term is likely to 
counteract much of this effect. So for all cases, /c ~ 2 is a reasonable value, but the variability in 
the density profiles and kinematics is sure to produce variation around this value. 



For our exploration of intrinsic galaxy quantities, we adopt a quantity, S = \/Kot + 2.00"^, 
which is just -v/2So.5 but with the benefit that it approaches Kot in the low-o" limit. This provides 
a quantity that can be directly compared with the velocity due to all mass, Vdrc- The constant k is 
a function of radius, and therefore it is useful to determine the range in radius for which A: ~ 2, or 
alternatively S ^ V^irc- In order to examine the value of as a function of radius, we stack all of the 
simulated galaxies, normalizing radii to R^q. For this analysis, we separate the particles into bins 
of width equal to 0.2i?5o. For each bin, we determine an average velocity and velocity dispersion 
for the stellar particles. This provides us with intrinsic quantities analogous to T4ot and a. These 
values are used to calculate iSintrinsic- V^irc = \/ GM{< r)/r is calculated for each bin by totaling the 
masses of all particles with a radius less than the bin, including baryons and dark matter. In Figure 
[6l we plot the ratio of S'intrinsic and Fcirc as a function of radius averaged over the velocity profiles of 
every snapshot within the simulation suite. One can see that at small radii there is no one-to-one 
correspondence between 5intrinsic and Vdrc as the ratio (5'intrinsic/V'circ) becomes closer to 2 or 3 on 
average with a large scatter between simulated profiles. However, at radii equal to or larger than 
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i?50) the ratio is quite close to unity and the scatter is about 30% in either direction. Therefore, at 
radii larger than R^q, ^intrinsic is a good indicator of total enclosed mass. For simplicity, we choose 
a radius of R^q at which to draw a value of ^intrinsic for comparison across the entire sample. 



5.2. Correlation between S'intrinsic and Vdrc 



Now that we have demonstrated that Sintrinsic is roughly comparable to Vdrc at R^q, it is 
illustrative to plot the relation between S'intrinsic and V^irc for the sample. We take these values at 
i?5o and plot them for 500 randomly chosen snapshots in Figure [71 If the two galaxies were separated 
by less than the sum of their stellar half-mass radii then they were treated as one. The x = y line 
represents equality between the two values. One can see that ^intrinsic is a good indicator of enclosed 
mass. The scatter in the relation increases with mass. This is the result of the distribution of galaxy 
types in the simulations. There are considerably more simulations run for the larger galaxies, and 
these simul ations have a wide va riety of orbits that can introduce significant spread in remnant 



properties (ICovington et al.l 120081 ). Therefore our simulation set will automatically produce more 
spread at higher masses. However, even for the larger mass cases, whose spread is likely more 
representative of the real universe than the low mass cases, log (S'intrinsic) ~ log(Vcirc) to typically 
within 0.15 dex. 



5.3. Correlation between ^observed and Intrinsic Vdrc 

A final question of interest is whether or not the observed quantity, S, can be used to estimate 
total galaxy mass, including the contribution from dark matter. To address this question, we plot 
the relation between ^observed and the intrinsic Vd^c in Figure[8l Again, we have moved the factor of 
2 in the So. 5 quantity over to the a in order to produce a correspondence between 5observed and Vdrc- 
Here we calculated a projected R50 using an elliptical aperture and measure a three-dimensional 
Vdrc within that radius. The V^t values are taken as the value of the fitted arctan curve at R^q 
rather than the asymptotic value. The observations are assumed to be at a redshift of one. As 
can be seen in Figure El there is a rough correspondence between ^observed and V^irc- Therefore, we 
take this result as a confirmation that 5observed can be used to track central galaxy mass, including 
contributions from dark matter. There is a tendency for ^observed to be a bit lower than Vdrc- 
This is likely the result of projection effects and counter-rotation, both of which would serve to 
reduce Sobservod compared with ^intrinsic- The outliers to low ^observed are typically overlapping 
cases where the kinematics are poorly represented by an arctan curve. The fact that the intrinsic 
relation is significantly better than the observed one suggests that it may be possible to improve 
the observational measurement of S. The assumptions of the arctan rotation curve and flat velocity 
dispersion often fail for the simulated galaxies, particularly when two galaxies overlap. Future work 
should examine possible improvements to the observational algorithm. 
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6. Discussion 



The observational results from iKassin et al.l (j2007l ) were interesting for several reasons. They 
demonstrated that a new kinematic quantity, So. 5, was capable of producing a relation that unified 
all galaxy types, including disturbed and merging cases. However, this study left open several 
questions that cannot be answered from the observations alone. First, what sort of mechanisms 
could produce the large scatter in the stellar mass TF while preserving the So.s-stellar mass relation? 
A second question is whether the relation is telling us something deep about the observed galaxies, 
or whether it is an artifact from the observations. Here we used a suite of galaxy merger simulations 
to address these questions. 

We find that galaxy mergers and interactions are capable of reproducing the observed relation. 
Specifically, there are merger stages during which the galaxies overlap and the rotation velocity 
significantly drops. This overlap results in increased dispersion which when included in a kinematic 
relation helps to reduce the scatter. Additionally, within a given merger one can see that rotational 
support is transformed into random pressure support. This transformation occurs in such a way 
that the scaling between log(So.5) and log(Mstar) is roughly conserved. Thus, galaxy mergers move 
galaxies up along the So.s-stellar mass relation, and the relation effectively unifies progenitors, 
mergers, and remnants. 

We examine the observational effect of blurring to see whether it could be responsible for 
the observed additional dispersion and lack of rotation. We find that the ROTCURVE algorithm 
performs well for the range of blurring found in the observations. We also find that S0.5 is typically 
conserved by blurring. Therefore it could be responsible for some of the reduction in scatter, but the 
effect is much too small to account for the small observed scatter in the So.s-stellar mass relation. 

Our examination of the kinematic scaling laws of the galaxies, and how they change with radius 
suggests that the appropriate constant in the S0.5 parameter is ~ 2. Additionally, we find that S 
at -R50 is a reasonably good tracer of V^irc at i?so- Therefore, it can be used to infer the galaxy 
mass within including contributions from dark matter. However, the details of this relation 
warrant more study. Specifically, the relation between S'lntrinsic and Vdrc is much better than that 
with S'observed, Suggesting that it may be possible to improve the observational measurement of 
S and consequently of the total enclosed mass. Additionally, it would be interesting to compare 
the merger si mulations with ava i lable integral held unit observations of nearby, z ~ 0.6, and 2; ~ 2 
galaxies (e.g. lArribas et al.ll2008l : lYang et al.ll2008l : iForster Schreiber et al.ll2006l . respectively). It is 
also likely that kinematics could be used as an indicator for merging galaxies, and possibly merger 
stage, since mergers often produce TF outliers and typically have significant asymmetry in their 
kinematics, but the development of new kinematic indicators of merging is left for future work. 
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Fig. 1.— 
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Fig. 1. — Time evolution of rotation velocity (red dotted line), velocity dispersion (blue dashed line), 
and So.5 (solid and long-dashed black lines) during a single merger simulation of two Sbc galaxies 
with stellar masses of ~ 5 x lO^^M© initially on a parabolic orbit. With each encounter, the rotation 
velocity decreases and the integrated velocity dispersion increases. However, So. 5 remains relatively 
constant throughout the merger, increasing primarily when the two galaxies coalesce. Care is taken 
to mock observe the simulations just as in the actual observations, assuming z = 1. The solid 
portion of the S0.5 line denotes the snapshots where only a single progenitor is observed in the slit, 
and the long-dashed portion denotes snapshots where both progenitors are obser ved in the slit. 



The merger simulation uses the smoothed-particle hydrodynamics code GADGET (jSpringel et al 



2001 



) and includes gas, star formation, and stellar feedback. The galaxy images are c alculated 
from the simulation using the dust radiative transfer code SUNRISE (jjonsson et al.ll2006l ) and are 
shown at four different times during the merger: 1) Before any encounters, 2) shortly after the 
initial encounter, 3) at the final coalescence, and 4) the remnant after the galaxies have coalesced. 
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Fig. 2.— 



- 23 - 



Fig. 2. — The TF relation (left) and So.s-stellar mass relation (right) for mock observations of 
simulated merging galaxies at redshifts of z ~ 0.03 (top), z ~ 0.3 (middle), and z ~ 1.0 (bottom). 
Each dot represents a mock observation of a single snapshot vie wed from a single angle. In order to 
compare with the observational results from lKassin et al.l (|2007l ). 500 such images were chosen from 
the simulation set at random. Red points are mock observations where both galaxies are present 
in the slit (i.e. close encounters and the merger remnant), and black points are mock observations 
of single simula t ed ga laxies. In the TF plots, the solid line is the high-reds hift TF ridge l ine fr om 



Conselice et al.l (|2005l ). with the dotted lines representing the scatter. As in lKassin et al.l (j2007l ). a 
significant number of the simulated galaxies scatter to low 14ot- In the observations, these galaxies 
have disordered or compact morphologies. Similarly, in the simulations, the majority of the cases 
scattered to low Vxoi are either undergoing or have recently undergone an encounter. In the Sq.s 
plots, the solid line is the fit to the observed So. 5 relation at z ~ 1.0, and the dotted lines depict 
scatter in the relation. The dashed line is the best fit to the simulations. Slopes (m), zero points 
(7/0)1 and scatters (a) for each fit are listed {y = yo + rnx). Including velocity dispersion greatly 
reduces the scatter and brings the progenitors, disturbed galaxies, merging galaxies, and merger 
remnants onto a single kinematic relation. The relation and scatter found for the simulated galaxies 
are comparable to the observed relation. Numbers on the plots show the location of the various 
numbered merger stages from Figure 1. The bottom panel has no '2' because it overlaps '1'. 
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Fig. 3. — Effect of blurring on kinematics. We take the four snapshots from our fiducial Sbc merger 
(undisturbed disk=top left, disturbed disk=top right, merger=bottom left, and remnant=bottom 
right) mock observed at z = 1, and analyze the kinematics after blurring the particle positions 
by gaussians with a's of 0.0", 0.1", 0.3", 0.6", 0.8", and 1.6". Blurring is normalized using the 
(unblurred) radius that contains 50% of the stellar mass (i?5o)- Lines and symbols are rotation 
velocity (blue triangles), a (red crosses), and So. 5 (black diamonds). Since the physical radii are 
different for each stage, and the blurring is for specific angular scales, each plot shows a different 
range of R^q/Rumv 
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Fig. 4. — The Stellar Mass TF and So.s-Stellar Mass relations analogous to Figure 2 with the 
analysis from the simulations done using only new star particles. 
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Fig. 5. — Histogram of the power law slope fit to the density profiles near R^q of each snapshot in 
the simulation set. 
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Fig. 6. — Ratio of 5intrinsic and Vdrc as a function of radius averaged over all simulation snapshots 
and normalized to R^q. Blue range shows the la scatter in the profiles. At radii larger than ~ R^q, 
Smtrinsic is a good estimator of total enclosed mass. 
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Fig. 7. — Relation between intrinsic S'intrinsic and Vdrc at R50 as observed in the galaxy merger 
simulations. Red points are cases for which both galaxies were analyzed together, most of which 
are remnants. The x = y line shows the rough equivalence of ^intrinsic and Vchc- 
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Fig. 8. — Relation between S'observed and Vdrc as measured at R^q. Red points are cases for 
which both galaxies are analyzed together. The line plotted is the x=y line, demonstrating the 
correspondence between the two values. 



